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Periodic Orbits on a Surface of Revolution. * 

By Daniel Buchanan. 



§ 1. Introduction. The object of this paper is to determine the periodic 
orbits described by a particle which moves, subject to gravity, on a smooth 
surface of revolution, the axis of which is vertical. Let us denote the equation 
of the surface by 

x 2 + y 2 -2pz + 2e<j>{x 2 + y 2 1 z)=0, (1) 

where p is a positive constant, e an arbitrary parameter, and <j> a power series 
in x 2 + y 2 and z, converging for x 2 + y 2 and | z | sufficiently small. There 
will be no loss of generality if we suppose that the constant term of <j> is zero, 
otherwise it could be eliminated by a linear substitution for z, that is, by a 
translation of the zy-plane along the z-axis. Since the equation (1) vanishes 
with x 2 + y 2 and z, it may be solved, by the theory of implicit functions, for 
x 2 -j- y 2 as a power series in z converging for \z\ sufficiently small. Hence 
the equation (1) may be expressed as 

F{x, y, z)=x 2 + y 2 — 2pz + 2ef(z) = 0, (2) 

where / is a power series in z converging f or | z | sufficiently small. 

For £= the equation (2) represents a vertical paraboloid of revolution. 
The generating parabola has its axis coinciding with the positive z-axis and 
its semi-latus-rectum equal to p. Periodic orbits are first constructed when 
the particle moves on this paraboloid. Then the analytic continuation of 
these orbits is made with respect to e and in this way the orbits for the more 
general surface (2) are determined. It was for this reason that the parameter 
e was introduced in (2). 

The form of the equation of the surface (2) was suggested by a somewhat 
similar equation used by Poincare in his memoir, Sur les Lignes Geodesiques 
des Surfaces Convexes: Trans. Am. Math. Soc, vol. vi (July, 1905). Geodesic 
lines play the same role in Poincare's memoir as periodic orbits do in this 
paper. 

It will be readily seen that when e = the problem here considered is 
somewhat analogous to the well-known problem of the spherical pendulum, f 

* Presented to the American Mathematical Society, Sept. 5, 1918. 

t For a complete discussion of the spherical pendulum, including the horizontal 
as well as the vertical motion, see Moulton's Periodic Orbits, chap. Ill, also Rendiconti 
del Circolo Matematico di Palermo, Vol. XXXII (1911), pp. 338-365. 
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§ 2. The Differential Equations. If the particle is of unit mass and 
moves without friction, then the differential equations of motion are 

x" = X, y" = Y, z" = Z — g, 

where X, Y , and Z are the normal reactions due to the surface, and g is the 
acceleration due to gravity. Since the surface is assumed to be smooth, the 
normal reactions at any point are proportional to the direction cosines of the 
normal at that point, and therefore the differential equations become 

x" = X = \F X = 2\x, -, 

y" = Y = \F y =2\y, t (3) 

z" = Z — g=\F ls — g = 2\{— p + £ / s )— g,\ 

where A is a factor of proportionality. 

These equations admit the vis viva integral 

x' 2 + y' 2 + z' 2 = 2g(c-z), (4) 

where c is the constant of integration. 

Since the equation of constraint (2) does not contain t explicitly, the 
factor A can be obtained by differentiating F( x, y, z) twice with respect to t 
and eliminating x", y", and z" from the result by means of the differential 
equations (3). Then on making use of the relations 

x * + y 2 = 2pz — 2ef, x' 2 -\-y' 2 = — z' 2 + 2g(c — z), 

from (2) and (4), respectively, we find 

2g(z-c)-g(p- eft) + ^ (1 - «f„) 
w ' ' f + 2pz-2e(f + pf g )+e 2 f/ ' l ' 



The part of 2X which is independent of e is 
= 2g(z — c)— gp + z n 



„/2 



(6) 



(A). Periodic Orbits on a Paraboloid of Eevolution. 

§3. Proof of Existence of Periodic Orbits. For e = the surface (2) 
becomes 

x 2 + y 2 — 2pz = 0, (7) 

and the differential equations (3) become 

x" = \ (s x, y"=X. y, z" = — p\. — g, (8) 
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where A is defined in (6). The last equation of (8) is independent of the 
first two and admits the integral 

„ 4gz(c — z)+c 1 
Z ' P + 2z ' (9), 

where c t is the constant of integration. 

If the vertical motion is to be periodic, z cannot increase indefinitely and 
therefore z' must vanish for some value of t, t = t , say. Suppose z = z at 
t = t . Then since z' = at t= t , it follows that the constant of integra- 
tion in (9) has the value 

c 1 = — 4 : gz (c — z ). 
Hence the integral (9) becomes 



4 



iff (z -z)(z -c + z). (10) 



p + 2z 



It is readily seen from (10) that z' vanishes for z = z , and also for 

2=C — Z . 

As p is assumed to be positive, no part of the paraboloid (7) will lie 
below the a;«/-plane. Then z and z cannot be negative and therefore, for real 
initial conditions, c cannot be less than z . Three cases arise according to the 
values assigned to c. They are : 

Case I. « = c < 2z . 



Case II. c=2z 
Case III. c > 2z, 



0- 



Case I. Suppose z ^ c < 2z or c — z < z . Then z' is real so long 
as z lies in the interval z =gl z jg; c — z . If the particle is started in the plane 
z = z 0> it cannot remain in this plane for z = z =?^= does not satisfy the 
z-equation of (8). Moreover, it cannot move above this plane since z < z , 
hence it must move below this plane. Consequently z decreases or z' is nega- 
tive, and the negative sign must be taken in (10). The particle continues to 
fall until it reaches the plane z = c ■ — z where the velocity again vanishes. 
As z = c — z is not a solution of the z-equation in (8) and as z cannot be 
less than c — z it must increase, or the positive sign must be taken in (10). 
The particle then rises until the plane z = z is reached where the velocity 
again vanishes and the radical in (10) changes sign. Hence the particle 
oscillates between the two horizontal planes z = z and z=c — z , the latter 
plane being the lower plane. 

4 
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Case II. If c = 2z or c — z = z , then 2' would be real only when 
z = z and would therefore be identically zero. In this case the particle 
would revolve in the horizontal plane z = z 0) the centrifugal force and the 
reaction of the surface being just sufficient to overcome gravity. 

Case III. If c > 2z or c — 2 > z it can be shown by an argument 
similar to that used in discussing Case I that the particle oscillates between 
the same two planes z = z and z = c — z , as in Case I, but the latter plane 
is the higher plane. 

Thus if c is positive and not less than z 0) the vertical motion of the 
particle is periodic, and therefore the last equation of (8) must admit a 
periodic solution. 

It could also be shown by the usual analytic existence proof that the 
last equation of (8) admits a periodic solution. It is not necessary, however, 
to establish the existence of a periodic solution by either method for, by 
Macmillan's theorem, quoted below, it is shown that if the fqrmal construc- 
tion of a periodic solution can be made, then this solution will converge under 
suitable conditions. Macmillan's theorem is as follows : * 

If 

-TT = fi(vi, • • • , x„; ii,t) (t = l, • • • , n) (11) 

is 'a system of differential equations in which the ft are expansible as power 
series in x u ■ ■ • , x n , and p., vanishing for x t = - • • = x„ = p. = 0, with coeffi- 
cients which are uniform, continuous, and periodic functions of t with the 
period 2jit: and if the ft converge for :£ t :£ 2 jit, when \ X) | < pj, \ p. | < r, 
then the solutions xi(t) are expansible as power series in p., or any fractional 
power of ji, which converge for all t in the interval ^ t ^ 2jt provided \ p. | 
is sufficiently small. If the constants of integration can be determined at 
each step so as to make the solution formally periodic with the period 2jit, then 
the solutions so determined will be periodic and converge for all finite values 
of t provided | /* | is sufficiently small. 

§ 4. The Vertical Motion in terms of Elliptic Integrals. The integral 
(10) may be written 



fV P + 2zdz 2V ^p =2V ^_g, (100 

We desire to express this integral in terms of elliptic integrals. The neces- 
* Trans. Am. Math. 800., Vol. XIII, No. 2, pp. 146-158. 
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sary substitution is found to be 

y/mn(l — kw) 

p +2z = =—7 -> 

1 + KW 

where 

m = p + 22 , « = p -\- 2c — 2z , k == — — —. 

Since p, c, and 2 are all positive, then m is also positive. If 2 < p/2 + c 
then m is likewise positive, and therefore k is real and less than 1 numerically. 
When the above transformation is made in (10') we obtain 

2y/mn f (1 — kw) dw 



C (1— «) ** n w~ (f f } 



Vm+ V«*' (1 — km>) V(l — w2 ) (1 — k2 w 2 ) 
which simplifies to 

2's/mn 
Vm -f- V« 
where 



{**-"'-tt^h£\-*w-u> 



/* dw 

1 ~ J V(l — » 2 )(1 — K 2 » 2 ) 

an elliptic integral of the first kind, and 



/,= 



J (1— K 



(1-A 2 ) V(l — w 2 )(l — * 2 w 2 ) 

an elliptic integral of the third kind. Thus t can be found for any given 
value of w or z. For most purposes, however, it is more convenient to have 2 
expressed as a function of t. This will be discussed in the next section. 

The half-period, T/2, of a complete oscillation can be obtained by 
integrating (10') between the limits z and c — z . The corresponding limits 
for w are — 1 and + 1. Hence 

2V 'mn f r+ % dw 






± V(l— W 2 )(l — K 2 W 2 ) 

/*+* dw 

_ J _ x (1 — K 2 W 2 ) V(l — W 2 )(l — K 2 W 2 ) 

2k r I 1 — w 2 ~i +1 ) 

1 — K 2 L\l — K 2 W 2 J. 1 J' 
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1+ 1 



since Wn^l*- ' then 

„ 4V«w f /•* dw 

" V7( VW + V«) K o V(l — w 2 )(l — «*w 2 ) 

J 1 Jw I 

(1 — K 2 W 2 ) V(l — W 2 )(l — K 2 W 2 ) )' 

§ 5. Construction of a Periodic Solution for the Vertical Motion. Let 

z = i; + c/2, (12) 

where v is the new dependent variable. Then the last equation of (8) becomes 

±gv + </ 2 

i;" = , 

Jc + 2v 

where fc = p + c. This equation can be expanded in the form 

v" = — {1/1) (4:gv + v' 2 ) [1 — 2v/h + 4tt 2 /& 2 

+ • • • (-2v/ky+ • • • ], (13) 

which will converge f or | » | < | -fc/2 |, or 2 < p/2 -j- c. Hence the particle 
must not rise above the plane z = p/2 + c, being the plane at the distance c 
above the horizontal plane passing through the focus of the generating para- 
bola of (7). 
Now let 



v _ yW ,t — 1 = IVI/J (1 + S) r, (14) 

where y and 8 are arbitrary parameters. When (14) is substituted in (13), 
the factor y can be cancelled off in the resulting equation, and we obtain 

w [(l + 8)w + ^ 2 ] [1 — 2p.w + - • •(— 2pu)l + - • • ], (15) 

where /t = y/Tc, and the dots denote derivation with respect to t. 

The terms of (15) which are independent of /a and 8 are w -\- w — 0, 
and the period of the solution is 2ir in t. Now an existence proof, based on 
Poincare's extension to Cauchy's theorem, would show that a periodic solution 
of (15) exists having the form, 

00 00 

w = % wax,*, 8 = 2 S,y , ( 16 ) 

where the Wj are periodic functions of t with the period %ir, and the 8,- are 
constants; and that this solution converges for all values of t in the interval 
^ t :S 2ir, provided | /«• | is sufficiently small. Instead of making this rather 
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simple existence proof, we shall assume the form of solution (16), make the 
formal construction, and then show that Macmillan's theorem will apply to 
prove the convergence of the solution. 
Since 

z = v + c/2 = yw + c/2 = hixw + c/2, (17) 

and since z' = at t = t , then at t = (or t = t ) we may choose, 

w(0) =a, w(0) =0, 

where a is a real constant. When these initial conditions are imposed on 
(16), we have 

w Q {Q)=a, w y (0)=0, i=l,- • -oo, w,-(0)=0, ; = 0,1,- • ■ oo. (18) 

Let (16) be substituted in (15) and let the resulting equation be 
denoted by (15'). This equation is an identity in /a and we may equate the 
coefficients of the same powers of /*, thus obtaining differential equations which 
define the various Wj. These equations are to be integrated, and the constants 
of integration and the various 8j are to be chosen so that the solutions will be 
periodic and satisfy the initial conditions (18). 

Step : Coefficient of p.". The terms of (15') that are independent of 
/* give the equation w + w = 0, and its solution which satisfies (18) is 
w = a cos t. 

Step 1: Coefficient of /*. The differential equation obtained from the 
terms in /a is 

w x -f w x = — 8^0+ 2w 2 — w 2 = — aS x cos t + (a?/2)(l + 3 cos 2t) .(19) 

In order that w^ shall be periodic the coefficient of cos t in the right member 
of (19) must be zero, otherwise the solution would contain the non-periodic 
or Poisson term rsinr. On putting 8 1 = 0, the general solution of (19) 
becomes 

w x = A x cos t + B ± sin t + (« 2 /2) (1 — cos 2t) . 

From the initial conditions (18) it follows that the constants of integration 
A x and B 1 are both zero, and the desired solution at this step becomes 

S t = 0, w x = (a?/2) (1 — cos 2t) . 

Step 2 : Coefficient of /i 2 . The differential equation at this step is 

w 2 -(- w 2 = — 8 2 w -\- 4w w 1 — 4w 8 -)- 2w v w 2 — 2w w 1 
= — ft(8 2 + y 2 a?) cos t — %a 3 cos 3t, 
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and the desired solution is 

8 2 = — %a 2 , w 2 = % 6 fr 3 ( — cos t + cos 3t) . 

The remaining steps of the integration are entirely similar to the pre- 
ceding step and by an induction to the general term it can readily be shown 
that the process can be carried on indefinitely. It will be observed that 8} is 
zero if j is odd, and that each Wj carries the factor a 1 * 1 and is a sum of cosines 
of multiples of t having the opposite parity of /. The highest multiple of r 
in Wj is O'+l). 

So far as the computation has been made, the solution of (15) is 

W = a COS T + (fM 2 /2) (1 — cos 2t) ^r- (cos t — cos 3t) + • • • , I 

8 iaY + ---. J 

The convergence of this solution will now be considered. In order to 
put the differential equation (15) in the same form as (11) of Macmillan's 
theorem, we let 

w = W v dWJdr = W 2 , 
then (15) becomes 

dW 1 /dr = f 1 (W 1 ,W 2 ;p,T)=W 2 , 

dW 2 /dr = f 2 (W 1 , W 2 ;p,r) = — ;[(1 + 8)W X + W 2 2 ][ l — 2 lt ,W 1 + • • • ]. 

The ^ and f 2 are obviously expansible as power series in W lt W 2 , and /*, and 
vanish for W t = W 2 = p = 0. The coefficients are constants and therefore 
satisfy the condition that they shall be uniform, continuous, and periodic 
functions of r. Further, f ± and f 2 converge for W u W 2 , and /t sufficiently 
small numerically. Hence equations (15) satisfy the conditions of Macmil- 
lan's theorem and consequently the solution (20) converges for | /* | suffi- 
ciently small. 

When (20) is substituted for w in (17) we obtain 

z = c /2 + fc|>acosr+^(l — cos 2t)—^ (cost — cos 3t)+ • • • J. 

As n and a are both arbitrary and as they occur only in products as indicated, 
we may suppress either without loss of generality. Let us suppose a = 1. 
Then the periodic solution for the vertical motion of the particle becomes 

z = e/2 + *D» cos t + (,i 2 /2) (1 — cos 2t) — ( V/16) (cos t — cos 8r) 

+•••], (21) 

8 1^ 2 +( )^ + - ... 
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In the expansion for z, the multiples of t in the coefficient of p.* have the same 
parity as / and therefore 

Z (ir) = c/2 — hp = c/2 — y. 
Further 

z(0) = c/2 + Tcp. = c/2 + y. 

Hence the particle oscillates between the planes z = c/2 + y and z = c/2 — y. 
The parameter y is therefore a scale factor denoting the amplitude of the 
oscillation on either side of the plane z = c/2. 

§6. The Horizontal Motion. When the solution (21) has been substi- 
tuted in the first two equations of (8) and the transformation 



t — to = W(Jc/g)(l + 8)r, 8 = -^+..., 
has been made, we obtain 

* + */P [i + «i/» + *•/** + ■ • - + M + - • •]*— 0, (32) 

and the same equation in y, where 

0i = — cos t, 2 = — % (1 — cos 2t), 

and the remaining 6j are likewise sums of cosines of multiples of t having 
the same parity as j, the highest multiple being j. 

Equation (22) is similar to the one first discussed by Hill * in 1871' 
in his celebrated memoir on the motion of the lunar perigee. A very complete 
list of references to the literature of the differential equations of this type is 
given by Baker on page 134 of his memoir " On Certain Linear Differential 
Equations of Astronomical Interest." f 

Three cases arise in the solution of (22) depending upon the values of 
(k/4=p).% They are 

Case I. i*Vk/p ^= and y/k/p not an integer. 
Case II. | \lk/p ^ and y/k/p an integer. 
Case III. \\lk~Jp~= 0. 

Since k = p -j- c and since p and c are both positive in the physical 
problem under consideration, it follows that \\Jk/p =^= and Case III need 
not be considered. 



*The Collected Works of <?. W. Hill, Vol. I, pp. 243-270; Acta Mathematica, 
Vol. VIII, pp. 1-36. 

t Phil. Trans, of the Royal Society of London, Series A, Vol. 216, pp. 129-186. 
t Compare also Moulton's Periodic Orbits, Chap. Ill, § 52. 
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Case I. — This may be regarded as the general case. 

The form of the solution of (22), first given by Floquet,* is 

x==e ia.T U> (23) 

where a and u are power series in n, the former having constant coefficients, 
the latter periodic coefficients with the period 2v in r. When (23) is substi- 
tuted in (22) we obtain 

■ii + 2tati — a 2 u + h/p[i + V+ V + • " * + 0>y + • • ' ]« = 0. (24) 

Now let 

« = u o + Uit* + w z/* 2 + ' • * > 

be substituted in (24). We obtain a differential equation (24'), which is an 
identity in /*. Since the right number is zero, then the coefficient of each 
power of /* must also be zero. On solving the various differential equations 
thus obtained, we determine the % and aj, the % as periodic functions of t 
having the period 2w, and the a, as constants so determined that the Uj shall 
be periodic. 

Since the solution (23) is later multiplied by an arbitrary constant, see 
equations (39), we may choose w(0) = 1, from 1 which it follows that 

w (0) = 1, Uj (0) = 0, ;= 1, • • • oo. (26) 

Step 0: Coefficient of p . From the terms of (24') that are inde- 
pendent of jx, we have the differential equation 



u + Wh/p w =0, 
and its solution is 

u = a + & e-*V*/* r } (27) 



where a and b are the constants of integration. Since W? is not an 
integer in this case, the term e~ iy/k/PT does not have the period 2ir and we put 
6 =0. Prom (26) it follows that a =l and the desired solution at this 
step becomes u = 1. 

Step 1 : Coefficient of (i. On equating to zero the terms in /i in (24') 

we obtain 

«i + iVk/p u\ = y/h/p di + h/p cos r. (28) 



* Annates de V&cole Normale Bupirieure, 1883-4. 
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The complementary function is 

k 1 =% + M'* t » 

«! and &j being arbitrary constants. 

It is well known in the theory of differential equations that the presence 
of any term in the right member which has exactly the same period as any 
term of the complementary function will yield Poisson terms in the particular 
integral, that is, terms containing t outside of trigonometric or exponential 
symbols. Since the constant part of the right member of (28) has exactly 
the same period as the constant in the complementary function, then non- 
periodic terms will arise in the particular integral unless a ± = 0. On putting 
Ql = the complete solution of (28) becomes 

-fc 

u 1 = a 1 + & ie -«Vfc/p t -|_ _ [cost — Wk/psinT]. (29) 

fC —— p 

From the periodicity and initial conditions it follows that 

p — fc 

and the desired solution at this step becomes 



[1 — cos t + Wh/p sin t], % = 0. (30) 



1 p — k' 
Step 2 : Coefficient of /* 2 . The differential equation at this step is 



k k 



u 2 +iVk/pu 2 = a 2 y/k/p + —- 2p{p _ h) (31) 

-j -, p- (cos t — Jcos 2t + i fyk/p sin 2t) — (k/p) cos 2t , 

and the solution which satisfies the periodicity and initial conditions is 
7pk 3 — 22p 2 fc 2 + 4p 3 fc — & 4 h* , • rrr • \ 

M2 = \p{p-kyl-±p) Tp-=W (cos T " * V h/p sm T) 



2pk 2 — 4p 2 k — k 3 . . . /TT - k(k + 2p) . . 

+ t-t T\-n rr cos 2r + 1 V k/p —. — 7N /T , x sm 2t, 

4,p(p — k)(k — 4p) 4(p — k)(k — 4p) 

02 = y 8 V^(f3)- 



(32) 



So far as the computation has been carried out it is observed that the 
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aj are real constants and that the Uj consist of sums of cosines and i times 
sines of multiples of r, the highest multiple being ;'. It will now be shown 
by induction that these properties hold in general. 

Let us suppose that a 1} • • • , an.!, u 0J u x , ■ • • , u n -x have been com- 
puted and that the various aj are real and that the Uj, j = 1, ■ • • , n — 1, 
have the form 

u s =k (A<-» cos h + iB<" sin It), (33) 

1=0 

where the A[ l) and B (l) are real constants. We wish to show that a n 

3 3 

can be determined as a real constant and that u n has the same form as (33) 

when / = n. 

Step n: Coefficient of p, n . The differential equation at this step is 



Un -\-Wk/pU„==.a n "VTc/p + Un(a 1} ■ • • , an-i.; U r , • ■ • , ttn-J. (34) 

The only undetermined constant which enters (34) is a» and it is written 
explicitly in so far as it occurs. The function U n is linear in u 1} • • • , u n _±. 
The' terms of U„ which arise from 2iaU in (24) have the form 

2iaiu n .i, I = 1, • • • , n — 1, 

and are cosines and i times sines of multiples of t, the highest multiple being 
n — 1. The terms which arise from a 2 u in (24) have the form 

— ( )aia m u„_ i i +m ->, l,m=0, ■ ■ ■ , n — 1, l-\-m-S^n — 1; 

where ( ) denotes 1 if I = m and 2 if l^=m. These terms have the same 
form as the preceding terms. When I = m = 0, the term — a 2 u n will cancel 
off with the term {k/^p)u n arising from the last factor of (24) and therefore 
the highest multiple of r in the sines and cosines is n — 1. The last factor of 
(24) gives, in addition to the term {Tc/^p)u n just considered, terms of 
the type 

(h/p) 6iUn-i, 1=1, ■ ■ • , n — 1. 

These terms are also of the same form as (33) but the highest multiple of r 
which they yield is n. Therefore 

U n = 5 [a (M) cos It + #><»> sin h], (35) 

1=0 l l 

where « <M) and & (n> are real constants. 

Let us now determine the periodic solution of (34). In order that the 
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solution shall be periodic, the right member of (34) must contain no constant 
terms. Hence 

o„ = — y/J/lca™, 
a real constant. The complete solution of (34) then becomes 

u n = a n + b n e- i ^T + U n , 

where a n and 6« are the constants of integration, and the particular integral 
U n has the same form as (35), except that it has no constant term. From 
the periodicity and the initial conditions we have 

&«=0, a n = — U„(0), 

and the desired solution for u n is the same as when j = n. Hence the prop- 
erties of aj and uj already stated hold in general. 

A second solution of (22) could be constructed in an entirely similar 
way but it is not necessary to make the construction as this solution can be 
obtained directly from the former solution. 

Let the solution of (24) already obtained be denoted by u(t, -\-i). 
Then one solution of (22) is 

x = e+ l "T u (t, +i). (36) 

Since the differential equation (22) is independent of i, a change in the sign 
of i in (36) will still give a solution, viz., 

x==e -ia Tu ( T) — i). (37) 

Thus a second solution can be obtained by changing the sign of i in the first 
solution. 

The determinant formed by the two solutions (36) and (37) together 
with their derivatives with respect to t is a constant,* and its value can be 
computed most readily when r = 0. This determinant is 



D 2[ia + m(0, + {)] = — i-\flc/p + terms in p, (38) 

which is different from zero for /* = and therefore remains different from 
zero for | /* | sufficiently small. The two solutions (36) and (37) therefore 
constitute a fundamental set and the most general solution of (22), as well 
as of the similar equation in y, is 

* Moulton's Periodic Orbits, § 18. 
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x = Ai_e iaT u w + A 2 e- iaT u i2 '> , 

5fc — p 



a 



-^[H^^-l 



i (1) = w(t, +t) = 1 -4- =r [1 — cos t + iy/k/p sin t]/* 

. r7pF — 22p 2 P + 4» 3 fc— -fc 4 F , . .-;—- . . 

, 2pk 2 — 4:p 2 k — k 3 

+ 7/ rrr; r^r cos * T 

4=p(p — k)(lc — 4p) 

k(k + 2p) 



,(39) 



+ <V*/P 4(p _ fc)( /_ 4y) ™ *)]** + ■ • • , 

where J. 1} J. 2 , J5 15 and I? 2 are arbitrary constants. 

By using Macmillan's theorem, quoted in § 3, it can readily be shown 
that the above solutions converge for all r in the interval ^ r iS 2ir provided 
| /* | is sufficiently small. 

Case II. Suppose y/k/p = v, an integer. Since k = p + c and p and 
c are both positive, then v > 1. Two sub-cases arise depending upon the 
value of v. They are : 

Sub-case I. v = 2, 
Sub-case II. v=^=2. 

Sub-case I. The construction proceeds as in Case I until equation (27) 
is reached, and this becomes 

Wo = flo + Ke- 2iT , (40) 

where both terms have the period 2ir. Prom the initial condition u a (0) = 1 
it follows that b = 1 — a , and the solution for u becomes 

u = a + (1 — ajr***, (41) 

with the constant a remaining arbitrary at this step. 

Step 1 : Coefficient of /*. When we equate to zero the coefficient of /* 
in (24"), that is in equation (24) after (25) has been substituted and \Jk/p 
replaced by 2, we obtain 

% + 2m = 2a 1 a — 2^ (1 — a ) <r 2iT + 2a e lT 

+ 2e-^ + 2(l — a )e- 3i T. (42) 
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The complementary function of this equation is 

«i = «i + b 1 e- 2iT . 

Since terms of exactly the same period as those of the complementary func- 
tion occur in the right member of (42), viz., constants and terms in e~ 2iT , 
the particular integral will contain non-periodic terms unless the constants 
and the coefficient of e~ 2iT are put equal to zero. Since we seek a periodic 
solution we therefore put 

2ajte =0, 2oj. (1 — Oo) = 0. (43) 

These equations are satisfied by a ± = 0, a arbitrary, but a must be different 
from 1 or this case would be the same as Case I. The complete solution of 
(42) then becomes 

«, = «, + & 1<r 2 *r _ %a e*r _j_ % e -ir _ % (i _ a<) ) e -u T . (44) 

Since «i(0) =0, then 

61 (%+tao). 

At this step a and % still remain undetermined. 

Step 2 : Coefficient of ju 2 . The terms of (24") which contain the factor 
H 2 give the equation 

u 2 + 2iu 2 = [2a 2 a + %a + 2] + [— 2a, (1 — a ) — %»o + 1 %]e" 2iT 

— 8(«i + %K 31 ' - *% (1 - '«o)e-* fT - 

On equating to zero the constants and the coefficient of e~ 2iT in the right 
member, as at the preceding step, we have 

2a 2 a — %to„ + 2 = 0, "1 

-2a 2 (l-a )-%a o + i% = 0.| ^ ; 

These equations are satisfied by 

„ 2 =±i/ 12 V217, a = 1 /2 + 1 /i4V317, (46) 

where the upper signs are to be taken together, also the lower. When equa- 
tions (46) are satisfied, the general solution for u 2 will be periodic, being 

u 2 = a 2 + b 2 e- 2i r — %a ± e { r + 5 / 12 « e 2 * r — %Go<r ir 
+ % (tei + % ) <r 3< ' + 5 /i2 (1 - o») e- 4iT - 

Since w 2 (0) = then 
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h= «2 + %«o — 4 %6- 

The undetermined constants at this step are a 2 and a t . 

Step 3 : Coefficient of i*?. The differential equation at this step is 

ii 3 + Ziu 3 = U 3 = — 2ia 3 u — 2ia 2 u\ -f- 2o 3 m + 2a 2 w x 

— u 6 3 + M i (4 — 4 cos % r ) + 4« 2 e ° s T - 

If we consider only the part of U 3 which involves the arbitrary constants that 
are determined at this step we have 

U 3 = 2a 3 a + a ± (2a 2 + %) + d x 

+ i-2(l-a )a 3 + ai (2a 2 -%)+a 2 e-^] (47) 

+ terms in e* 1 *, e 2 ^ e *3i T) e ±s* T) e -« T) 

where ^ and d 2 are known constants. To make u 3 periodic we equate to zero 
the constants and the coefficient of e~ 2ir in (47) giving 

2a 3 a + on (2a 2 + % ) + d t — 0, "1 

— 2(1 — a,)a 3 + ai (2a 2 — %)+d 2 — O.J <■ ' 

The determinant of the coefficients of a 3 and a^ in (48) is 

D 1 =±%V217, (49) 

where the upper sign here is to be taken with the upper signs in (46), simi- 
larly for the lower signs. Since this determinant B 1 is different from zero, 
equations (48) can be solved for o 3 and a t> and the solutions are unique for 
either set of values in (46). When (48) is satisfied, the solution for u 3 will 
therefore be periodic. Two constants of integration, a 3 and o 3> say, will 
arise from the complementary function. One constant, b 3 , say, will be deter- 
mined from the initial condition w 3 (0) =0, while a 2 and a 3 remain undeter- 
mined at this step. 

The succeeding steps of the integration are entirely similar to step 3 
just considered. So far as the constants of integration are concerned one of 
them, I) say, is determined by the initial condition Uj(0) =0 at the step 
where it arises. The other constant a, of the step /, is not determined until 
the step / + 2 is reached, where, on equating to zero the constants and the 
coefficient of e~ 2iT in the right member of the differential equation in M/ +2 , 
two linear equations in a; +2 and a } - are obtained and the determinant of their 
coefficients is the same as D t in (49). Hence <y +2 and a s can be uniquely 
determined provided one set of values is taken in (46). 
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It is thus evident that two solutions of (24") can be obtained in the 
foregoing construction, according as the upper or lower signs are taken in 
(46). The corresponding solutions of (22) are 

x = Xl = 8«C»ViiVaf ^+ . . . )r [ (i/ 2 _ i/ 14 V317) 

+ (% + Vi4 V%ri)e-^ + ( )/»+••■], 
and 

x = z 2 = e «(i- 1 A 2 vH?^+ . . • )r [(% + i/ 14 V217) 
+ (y 2 -yi4V217)e- 2 ^+ ( )/* + •••]. 

If we proceed as in Case I to obtain two solutions by changing the sign 
of i in the preceding solutions it would appear that two additional solutions 
could be obtained, but this is impossible since the differential equation (22) 
is only of the second order and admits of but two solutions. This apparent 
difficulty is overcome if the factor e iT is multiplied into the part of the solu- 
tions contained in the square brackets [ ]. Thus 



Xl = e W^t M »+ . . . ) T [ ( y 2 _ % 4 V217) e^ + ' 

+ (% -M4 4 V2i7)«-*r+ ( )^-f— • ], 

x 2 = e -«(V u V2iT „'+... )T [-(i^ + y 14 V2i7) e *r _j_ 
+ ( 1 /2- 1 /i4V317)^+ ( )/* + •••]. 



(50) 



It is obvious, therefore, that the two solutions x x and x 2 differ only in the sign 
of i or of V217. 

If the exponentials in the square brackets of (50) are expressed in 
trigonometric form, these solutions become 



x t = e ia ru w , x 2 = e- iaT u™, 

q^y 12 V217^ + - • • , 

m< 1 > = cost — t%V217sinT+ ( )/* + 

u™(+i)=u w (—i). 



(51) 



The terms represented by u m and w <2) differ not only in i but also in the sign 
of V217. They are power series in /*. with sums of cosines and tV2lT times 
sines of multiples of t, the highest multiple in the coefficient of /*' being 

y+i. 

The determinant of the two solutions in (51) together with their deriv- 
atives is a constant, as in Case I, and its value at r = is 

P 2 = -2[to + # 1 >(0)]=t2/ 7 V217+( )/* + • • • , 
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which is different from zero for /*. = 0, and therefore remains different from 
zero for | /* | sufficiently small. Hence the two solutions (51) constitute a 
fundamental set, and the most general solutions of (22) and the similar 
equation in y are 

x =, A^aru™ + A 2 e- ia ru m , y = B x e ia Tu m + B 2 e- ia ^u (2) , (52) 

where a, u a) , and m (2) are defined in (51), and A 1} A 2 , B u B 2 are the constants 
of integration. 

Sub-ease II. When v is an integer different from 2, the construction 
is the same as in the preceding sub-case until equations (43) are reached. 
The equations analogous to (43) are 

va^ — 0, va^fl — a ) = 0. (53) 

These are satisfied by % = 0, a arbitrary, but a must be different from 1 as 
in the preceding sub-case. 

The equations which correspond to (45) are 



(54) 











«o{ 


. v 2 (5v 2 - 
" a2+ 8(v»- 


^H<1 








(1- 


«o) { — 


Va2+ 8(v 2 - 


-,V}- 


Since 


a y^l 


in 


Case II 


, then the solutions of 


(54) are 








ftn 




v(5v 2 — 1) 





8(^-1) " ( 55 > 

The equations analogous to (48) are 

<*i { ™ 2 + ^(v 2 -l) } = °' ~ ms + " 2 ^ 0> = °' (56) 

where ( 3 O) is the constant part of 3 in equation (22). The solutions of these 
equations are 

Oi — o, a 3 = v0< o) . 

The constants a 2 and a 4 are determined in the same way at the next 
step as a 1 and a 3 were found in the preceding step. They have the same 
coefficient in the equations analogous to (56) as a ± and o 3 have in (56), and 
similarly for the succeeding steps. It will be found that 

% = 0, j/'=0, 1, • • • , v— 3, 

and that the remaining a, are, in general, different from zero. 
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So far as the computation has been made we find 
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,,2 r~(,-Up-l>T (T<I»«'T-| 

Wz L(v 2 — l) 2 ^8(v 2 — 4) J 2(v 2 — 1)L v — 



i(y-l)T 

r 






■ + : 

% = 0, a 2 = 



-2 



v + 2 



J. 



v(5v 2 — 1) 



8(v 2 — 1 

When these terms are substituted in (25) and the result in (23) we obtain 
one solution of (22). Let it be denoted by 

x==e iar u w, (58) 

Another solution can be obtained by changing the sign of i in (58), thus 

x = e-^Tu™, (59) 

where m (2) (+ i) = u w ( — ■£) . 

The determinant of these two solutions and their first derivatives is 

— 2[ii + fi<"(0)]— *+( )/» + •••, 

which is different from zero for \p\ sufficiently small. The solutions (58) 
and (59) therefore constitute a fundamental set, and the general solutions 
of (22) and the corresponding equation in y are of the same form as (39) 
or (52), viz., 

x = 4 1 e iOT w (l) + 4 2 e-<<"- w (2) ( y = B x e i <"u m -\- B 2 e- ia ru m , (60) 

where a and m (1) have the values in (58). If the exponentials in (58) are 
expressed in trigonometric form then 

«<» = SJ [iw + £«> cos(v ± l)r + iC7«' sin(v ± Z)t]^, (61) 

where A w , B ( /', and 0'/' are real constants. In particular 

A<» = 0, j = 0, 1, 2, B«> 0«), ; = 0, 1, 2, J = 0, 1, 2, 

# <0) = 1, B w = 



J5W= 

- 1 2(v 



v 2 — 1' 

v 2 
, fid) = _ 

— 1)' - 1 



2(v+l)' 
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B ™. 


(v 2 


— I) 2 '8(i/ 2 - 


-4)' 


#w 


B i2) ' 

+1 




v 4 




n(2) 




3 ( V + 1)( V 2_ 


1)' 


-2 


R(2) . 




v 2 






^+2 


32 


:(v+2f 




Prom 


(57) 


it follows that 










A^ = 0, 


/ — 


.0,1, 



2(v-l)(v 2 -l) ' 



32(v— 2)' 



The remaining ^.°' ) are different from zero, in general. 

If the factors e i( " /2)T and e~* ( " /2)T are taken with w (1) and w (2) , respect- 
ively, in (60), then w (1) and w (2) will have cosines and sines of (v/% ±l)r 
instead of (v ± l)r as in (61). 

§ 7. The Arbitrary Constants. Bach of the solutions for x, y, and z 
contains two arbitrary constants, and the values of these constants for the 
physical problem will now be considered. 

In this discussion no discrimination will be made between the two cases 
of the preceding section. It was for this reason that the same notation was 
chosen for the solutions (39), (52), and (60). 

Besides satisfying the differential equations (8), the solutions for x, y, 
and 2 must also satisfy (7), viz., 

* 2 + f — 2f« = 0, (7) 

and their derivatives with respect to t must satisfy the vis viva integral (4) . 
When the transformation 



t-t = iy/h/g(l + S)r, 8 = -^ 2 + - • • 

is made in (4), this integral becomes 

x 2 + y* + z 2 =(k/2)(l + 8)(c — z). (62) 

Since any horizontal section of the surface (7) is a circle with centre 
on the 2-axis, the xy-axes may be rotated about the 2-axis without changing the 
form of (7), and we may therefore suppose that the particle lies in the a;2-plane 
at the initial time, or y = at r = 0. Since m (1) (0) = u (2) (0) = 1, then it 
follows that B x + B 2 == 0. 

As 2(0) =0, the particle must be initially projected in a horizontal plane. 
Since any horizontal section of the surface of constraint is a circle with centre 
on the 2-axis and since w(0) = 0, then it follows that £(0) = 0. When this 
condition is imposed on the solution for x, we have 
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[t. + t* (1> (0)] [A 1 — A 2 ] = 0. 

The factor [to + it (1) (0)] is different from zero for J ju. J sufficiently small and 
therefore 

A x = A 2 . 

Hence the solutions for x and y become 

x = A[e i «ru w + e-^Tu™], y = B[e i ^u w — e- iar u™~], (63) 

where the now unnecessary subscripts on the constants A and B have been 
dropped. 

On putting t = in (63) and (21) we get 



x(0)=2A, a;(0)=0, 

2/(0) =0, y(0)=2B[ia + 4^(0)1, 

2(0)=c/2 + fc/*, «(0)=0. 

When these values are substituted in (7) and (62) we find 
A = ± i Vp(c + 2fyx) — ± iVp(c + 2y), 



(64) 



B = ± 



V fc(l + S)( C -2y) 
4[ta + « (l >(0)] 



Vfc(i+a)( c -2 Y ) 

4[a — BWfO)] 



Thus 4 is real and B is purely imaginary. Hence the solutions of equations 
(8), in terms of t, are 



x=± ly/p{c + 2y) [e* OT « (1) + e- ic »-w< 2 >], 

„ = _ ^V^(l+g)(c — 2y) r e * T U <i) e -* 0T w (2) l 

9 "*■ 4[a — tfi<«(0)] 

2 = c/2 -f- fc[/A cos t + i/* 2 (l — cos 2t) 

— %6/* 3 (eosr — cos3r) +• • • ], 



(65) 



Tc = p -+- c. 

The terms a, w (1) , and u i2) are defined in (39), (52), and (60) according to 
the value of y/k/p. The double signs in (65) depend upon the octants of 
space into which the particle is initially projected. 

If the exponentials in (65) are expressed in trigonometric form, the 
solutions of (22) for x and y become 



X=y/p(c + 2y) [{ } lC osaT— { } 2 sinaT], 



y = 



Vfc(l + 8)(c — 2y) 
2 [a — tti^O)] 



[{ } 1 sinor+ { } 2 cosar], 
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where { } t and { } 2 are power series in /* which have different values for the 
different cases of § 6. For Case I, 

00 j 

{ } 1 = l + 2 S[>/jcosZt]/*', 

)=1 1=0 
00 j 

{ } 2 = 2 S [&** sin ZtV, 
i=i i=a 

where a$i and bji are the coefficients of cos h and sin Zt, respectively, in Uj, 
equations (30), (32), and (33). 
For Case II, 

00 j 

{ } x = 5 2 [4<'> + B«> cos(v ± 1)t]iJ, 

)=0 1=0 ±l 

{ } 2 = S i[O«>ain(v±0r]^. 

/=! 1=0 *« 

The constants A^', JB W , and C ( '> have the same values as in (61) for sub- 
case II, and for sub-case I they are, in so far as the computation has been 
carried out, 

i.< 0) = 0, B^ = 1, 0< 0) = — % V»17. 

There still remain in the solutions (65) two constants which have not 
been determined. They are c and y (=&/i). Since c is the constant of 
integration arising in the vis viva integral, its value depends upon the initial 
velocity of the particle. As x' = z' = at t = t , the initial velocity may be 
denoted by y' . Then from the vis viva integral (4) it follows that 

y>*=2g(c — z )=g(c — 2y), or c — 2y = y'*/g. 

Since, further, c + 2y = 2z , then 

c = *o + *C/2<7, y-i(*o — #/»?)• 

Thus the constants c and y are functions of the initial velocity and the initial 
height. 

Let us now return to the three cases of § 3. Since z = c/2 + y these 
cases become : 

Case I. 0<yg c/2, 

Case II. y = 0, 

Case III. y < 0, | y | ^ c/2. 

In Case I and Case III the particle oscillates between the horizontal planes 
z = c/2 + y and z = c/2 — y. The two orbits are geometrically the same 
but the motion in the one orbit is half a period ahead of the motion in the 
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other. If y = c/2, then the initial velocity is zero and as there is no lateral 
projection the particle will move in the vertical parabola x 2 = 2pz. If 
y «= — c/2, then z = 0, and the particle is projected from the lowest point 
with the initial velocity y' — \/2gc. It will therefore move in the vertical 
parabola y 2 — 2pz which is dynamically the same orbit as x 2 = 2pz. The 
orbits when y = ±c/2 correspond therefore to the simple pendulum. 

In Case II when y = 0, then z = c/2, a constant, and the differential 
equations of motion become 

*"+ (g/p)*-o, y"+(9/p)y = o, *"-o. 

Their periodic solutions which satisfy the initial conditions are 

x = \/pc cosy/ g/pt, y — y/pc siny/ g/pt, z — c/2. 

In this case the particle moves in a circle the plane of which is parallel to the 
a;y-plane and at a distance c/2 above it. 

In the first paragraph of § 5 it is stated that the inequality z < p/2 -\- c 
must hold in order that the expression in (13) will converge. Now the 
maximum value of z is 

*-c/8+|y|, 

and therefore the above inequality becomes 

c/2 + \y\<p/2 + c, 
or |y|<*(P + c). 

Since | y | must not exceed c/2 in order that the initial velocity shall be real, 
it follows that the inequality z < p/2 -J- c will always be satisfied for real 
initial conditions. 

(B). Periodic Orbits on a Surface of Revolution. 

§ 8. The Method of Solution. — Let us next consider the construction 
of the periodic orbits described on the more general surface of revolution 
represented by the equation 

F{x,y,z)=x 2 -\-y 2 — 2pz+2ef(z) = 0. (2) 

The differential equations of motion have already been found, equations 
(3) and (5). The method of constructing the periodic solution of these 
equations is to first make the analytic continuation with respect to e of the 
periodic solution for the vertical motion obtained in (A) where « = 0, and 
then substitute this solution for z in the first two equations in (3). We thus 
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obtain two differential equations having periodic coefficients somewhat similar 
to (22). 

§ 9. The Equation of Variation. — Let us substitute in the last equation 
of (3) ~ 

8—*. + ..., i < M > 

where z denotes the solution obtained in (21), and £ is a function of t which 
vanishes with e. We obtain 

'i + ®J + ®^ = Z + e Z 1 + - •■+,%• + • • • , (67) 

where the undefined terms have the following properties : 

(1). The functions © x and © 2 are periodic functions of t, but we shall 
show that it is not necessary to know the explicit values of these functions in 
order to solve (67). 

(2). The functions Z , Z x , • • ■ , Z } denote power series in £ having 
coefficients which are power series in fx. with sums of cosines of multiples of t 
in their coefficients. These functions also contain additional terms in £ and | 2 , 
the former being multiplied by power series in ju. with sums of sines of mul- 
tiples of t in the coefficients, the latter by similar series except that they 
contain cosines. The function Z contains no linear terms in £ or £. 

If we neglect the right member of (67) we obtain 

|'+©J + © 2 |=0, (68) 

which is called the equation of variation. The generating solution is 2 = 3", 
or the expression for z in (21). 

Now it has been shown by Poincare * that if the generating solution 
contains an arbitrary constant which does not occur in the original differential 
equations of motion, viz., equations (3), then a solution of the equations of 
variation can be obtained by differentiating the generating solution with 
respect to this constant. Three constants occur in (21), viz., c, m, and /„, the 
latter entering implicitly through t. The constant c occurs in 2A, equation 
(5), and therefore enters the differential equations (3). The two remaining 
constants, n and t , do not occur in (3) and therefore may be used in applying 
Poincare's theorem. Each constant yields a solution of (68) and since the 
differential equation is only of the second order, both its solutions may be 
obtained from the generating solution. Hence it is not necessary to know 
the coefficients ® t and ® 2 in (68) in order to solve the differential equation. 

* Les MSthodes Nouvelles de la Mecanique Celeste, Vol. I, chap. IV. Loc. cit. 
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Consider, first, the constant t . Then one solution of (68) is 



«_*/»«,— 8 V^r+iT m/dt ^ 2 Vri+8) ^ (t) ' (69) 

where 

S(t) = sin T — ft. sin 2r—y iel i 2 (sin t — 3 sin 3t) + • • • . (70) 

Since this solution is later multiplied by an arbitrary constant, see equation 
(76), we may drop the constant factor of S(t) in (69) and take 

£ = £(r) (71) 

as the solution. 

Considering the constant /*, we obtain as the second solution of (68) 

where the parentheses ( ) denote differentiation in so far as /*. occurs explicitly 
in z. Now 

dz 

■x- =1cC(t) =&[oost + /x(1 — cosSt) 

— ~f (cos r — cos 3t) + • • • ], }■ (73) 

^ = — i jqjg, g|g = til— l + power series in /t 2 ]. 

Therefore the solution (72) becomes 

£ = fc[C(r)+TZS(r)], 
where 

IT = ii, 2 [ — l + power series in ja 2 ]. 

As in the previous solution we may drop the constant factor h and take 

£ = C(t)+tKS(t), (74) 

as the second solution of (68). 

The two solutions (71) and (74) constitute a fundamental set, since the 
determinant of these two solutions together with their derivatives is different 
from zero for | /* | sufficiently small, being 

A = — l + power series in /*. (75) 

Hence the general solution of (68) is 

t-N l 8(r)+N a lC(.r) + tEB(t)1, (76) 

where N t and N 2 are arbitrary constants. 
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From the way in which S(r) and C(r) were derived, it is readily seen 
that the coefficients of sin(; + 1)t and cos(j + 1)t in the coefficients of pi 
in S(t) and C(t), respectively, are equal. 

§ 10. Construction of the Solution for the Vertical Motion. 
Let 

oo 

l = S &, (77) 

i=i 

where the t) are to be periodic with the period 2w in t. 

Since the initial time was chosen for £ = so that 2(0) =0, it may 

now be chosen so that i:(0)=0 when e^O. Now 2 = 2 + £, and as 

1(0) = 0, it follows that |(0) = 0. Hence 

li(0)=0, j=l, ■ ■ ■ oo. (78) 

Let (77) be substituted in (67) and let the resulting equation be 
denoted by (6 7'). It is an identity in e, and on equating the coefficients of 
the same powers of e we obtain sets of differential equations which can be 
integrated and the various constants of integration can be chosen, as we shall 
show, so as to satisfy the periodicity and initial conditions. 

Two types of series occur in this integration and they are similar to 
S(r) and O(t) in (71) and (73), respectively. These series will be 
denoted by 



or 



B,(t), S^(r), S,(r), SW(t), j=l,2, 
C,(r), 0<»(t), C J (r),C^(r),j = l,2, 



according as they are similar to S(r) or C(r), respectively. 

Coefficients of e. When the coefficients of <■ to the first power are equated 
in (67') we obtain 

£ + ®A + ® 2 Ii = C (1) (t). (79) 

The complementary function of this equation is the same as the solution of 
the equation of variation, viz., 

^=n^S(r) +n»'[C(r) +rfS(r)], (80) 

where n™ and n™ are the arbitrary constants. By employing the method 
of the variation of parameters to find the complete solution, we have 

n^S(r)+n^[C(r) +tKS(t)-\ = 0, "| 

. . Y (81) 
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The determinant of the coefficients of n^ and hf in the above equations is 
the same as (75), and since it is different from zero equations (79) can be 
solved for w (1) and n™ . Thus 



2 



h™ = (1/A)0 (1 >(t) [C(t) +tKS(t)], 



(82) 



«a>= (1/A) C (1) (t) S(r). 

Since the coefficients of the same power of /* in (7 (1) (t) and C(t) are sums 
of cosines of the same multiples of t, the product C w (t)C(t) will yield, in 
addition to periodic terms, a constant 

wi = /*L>r +p™t* + • • ■ + PiV + • • • l 

where jp w are real constants. Then the integration of (82) gives 
n w = # <i) — |> T ft + £ (1) (t) + r Z <T (1) (r) ], 



(83) 
w a)=^(i)+<7a)( T ), 

where N™ and iV 2 (1) are the constants of integration. Since the coefficients 
of the sines and cosines of the highest multiples of t in the coefficients of /»' 
in $(t) and O(r), respectively, are equal, then it follows that 8 w (r) and 
(1) (t) have the same property. When (83) is substituted in (80) the 
complete solution of (79) then becomes 

£ 1= = N ™S(t) +N?IC(t) +tKS(t)1- ( it Pi S(t) + (l»^(r), 

where 1 (t) contains no terms independent of p. In order that £ x shall be 
periodic iV (1) must be given the value 

N (1) = — - = — (power series in ft), (84) 

and from the initial conditions (78) it follows that 

jy»» = 0. 

Hghog 

^x=(1»C 1 (t), (85) 

where 1 (r), like (Mt), contains no terms independent of /t. 

Coefficients of e 2 . Equating the coefficients of e 2 in (67) gives the 
differential equation 

L + ®Jz + ®A = (1»C (2) (r). (86) 

Except for the factor 1/ju. in the right member, this equation is similar to 
(79). The general solution is obtained in the same way as at the preceding 
step and is found to be 
6 
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S a -N™B(r) +N^[_G(t)+tKS(t)-\ +rp 2 S(r) + (1/ M )(J< 2 ' (r), 

where iV (2) and N<- 2) are the constants of integration, and p 2 is a power 
series in p. In order to satisfy the periodicity and initial conditions we 
must pnt 



2V< 2) =£?-= — (power series in /*), 2V< 2) =0. 

2 K p 2 ^ 1 

Hence 

£ 2 =(l/p 2 )C 2 (r), 

where 2 (t), contains no terms independent of p. 

The remaining steps of the integration can be carried on in the same 
way and by an induction to the general term it can be shown that 

&= (1>»)C.(t), 

where C n {r) as at the previous steps, contains no terms independent of p. 
On substituting the solutions for the various £y in (77) we obtain 

f-S C,(T)(e/rf/. (87) 

By applying Macmillan's theorem, quoted in § 3, it is found that the solution 
(87) converges for | c [ sufficiently small. 

When (87) is substituted in (66), the solution for the vertical motion 
becomes 

CO 

«-i+SO,(T)(«/(.)'. (88) 

Since e and p are both arbitrary we may put e = pp and therefore (88) becomes 

00 

3 = 2+2 C,-(tV. (89) 

This solution converges for | p | and | p | sufficiently small. 

§ 10. The Horizontal Motion. When the first two equations of (3) 
are transformed by the substitution 



t — *„ = |V(V3')(1 + <S) r, 8 = — &* + •••, 

already used in (A), and the value of z obtained in (89) is substituted in 2\, 
these differential equations become 



x + Oo + <?xP + 4> 2 p 2 + • • ■ ] « = 0, 

y + Oo + <l>ip + <p 2 p 2 + • • ■ ] y = o 



|' } 0M>) 
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where each $j is a power series in ju. with sums of cosines of multiples of r 
in the coefficients. The function <£ has the same value as the coefficient of 
x in (22). 

Now let p = a/x. Then the coefficients in (90) can be rearranged as 
power series in /* and the differential equations take the same form as (22). 
Hence the solutions are 

x = L ± e*/ 3 ^ -j- L 2 e^P^, y = M x e t P T v 1 + M 2 e-^ T v 2 , (91) 

where L 1} L 2 , M 1} and M 2 are the constants of integration, and ji, v t , and v 2 
are similar in form to o, u ±) and u 2 , respectively, of (39), (52), or (60) 
according to the value of yfk/p. 

When the initial values £(0)=y(0)=0 are imposed on (91), it 
follows that 

W + *i(0)][£i-£J-0, M. + M^O. 

Since i/3 + ^(0)^= for | p | sufficiently small, then 

L 1 = L 2 = L, say, 

and, from the second equation, 

M 1 = — M 2 = M, say. 

Hence the solutions (91) become 

x = Lle i ^v 1 + er^Vtl, y = M-[e*PT Vl — e"*/ 3 ^] . 

The constants L and M can be determined as in § 7, and it is found that L is 
real while M is purely imaginary. Suppose M = iN. Then 

00 j 00 j 

a; = 2L[cosy3T{l + 2 S/tWcosZr} — sin^r{S 5 /Wsin Zr}], 

.7=11=0 y=ll=l 

» j 00 j 

y = 2iV[sinMl + S S/t%i«»cosZT} — sin/? T {S S/tW'sin h}], 

where «<>'> and &<>> are real constants. These equations together with equation 
(89) are the solutions of the equations of motion of the problem under 
consideration. 
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